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Abstract 

In this paper wc try to find examples of integrable natural Hamil- 
tonian systems on the sphere with the symmetries of each Platonic 
polyhedra. Although some of these systems are known, their expres- 
sion is extremely complicated; we try here to find the simplest possible 
expressions for this kind of dynamical systems. Even in the simplest 
cases it is not easy to prove their intcgrability by direct computation 
of the first integrals, therefore, we make use of numerical methods to 
provide evidences of integrability; namely, by analyzing their Poincare 
sections (surface sections). In this way we find three systems with pla- 
tonic symmetries, one for each class of equivalent Platonic polyhedra: 
tetrahedral, exahedral-octahedral, dodecahedral-icosahedral, showing 
evidences of integrability. The proof of integrability and the construc- 
tion of the first integrals are left for further works. As an outline of 
the possible developments if the integrability of these systems will be 
proved, we show how to build from them new integrable systems in di- 
mension three and, from these, superintegrable systems in dimension 
four corresponding to superintegrable interactions among four points 
on a line, in analogy with the systems with dihedral symmetry treated 



1 



in a previous article. A common feature of these possibly integrable 
systems is, besides to the rich symmetry group on the configuration 
manifold, the partition of the latter into dynamically separated regions 
showing a simple structure of the potential in their interior. This ob- 
servation allows to conjecture integrability for a class of Hamiltonian 
systems in the Euclidean spaces. 

1 Integrable and superintegrable systems 

It is not easy to include all known completely integrable systems within a con- 
cise definition or even a list; they appear in almost every field of mathematics 
and mathematical physics, from algebraic geometry to supersjTiimetry theory 
touching classical and quantum mechanics, elliptic curves, minimal sufaces, 
number theory, Riemann surfaces etc [13], [3], [H], [23] [22] • In few words, by 
quoting F.Helein [13]: "...working on completely integrable systems is based 
on a contemplation of some very exceptional equations which hide a Platonic 
structure: although these equations do not look trivial a priori, we shall dis- 
cover that they are elementary, once we understand how they are encoded 
in the language of symplectic geometry. Lie groups and algebraic geometry. 
It will turn out that this contemplation is fruitful and lead to many re- 
sults" ; the word " Platonic" must be intended here in its purely philosophical 
meaning. In the following we consider classical time-independent completely 
integrable Hamiltonian systems. In this context, complete integrability co- 
incides with Liouville integrability: a Hamiltonian system with n degrees of 
freedom is Liouville integrable if it admits n functionally independent and 
Poisson-commuting first integrals (i.e. constants of the motion), including 
the Hamiltonian function itself; the integral curves of the system are then 
determined by quadratures |lj . Finite dimensional Hamiltonian systems can 
admit up to 2n — 1 functionally independent first integrals. Liouville inte- 
grable Hamiltonian systems with extra independent constants of the motion 
are said to be superintegrable [27]. A well known example of completely in- 
tegrable and superintegrable system is the Kepler system in the plane, whose 
Hamiltonian in polar coordinates (r, ?/') is H = ^{Pr + ^P^) + ^ where k is 
a real constant, its functionally independent first integrals are the Hamil- 
tonian, the angular momentum and one of the components of the Laplace 
vector (the other depends functionally on these three) [IT]. The integral 
curves of (time independent) Hamiltonian systems belong to the phase-space 
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(cotangent bundle, in the language of symplectic geometry), whose dimen- 
sion is twice the degrees of freedom of the Hamiltonian; the existence of 
2n — 1 first integrals implies that each integral curve is geometrically de- 
termined by the one-dimensional intersection of all the 2n — 1 dimensional 
hypersurfaces corresponding to the constants of the motions determined by 
the initial conditions, then, without need of integrating the Hamilton differ- 
ential equations, even if it can be practically impossible its explicitation in 
that way due to the mathematical complexity of the functions involved. If 
the independent constants of the motion are 2n — 2, then each integral curve 
will stay in their two-dimensional intersection and so on. The n Poisson- 
commuting independent first integrals Ki of a Liouville integrable Hamilto- 
nian, under the assumption of completeness and excluding eventual critical 
points of the Hamiltonian fields, generate the foliation of the phase space 
into the n-dimensional submanifolds determined by the intersection of the 
manifolds Ki = const., which are diffeomorphic to M'^ x for < k < n, 

where T"~'^ is the n — k dimensional torus (known as Liouville torus). If the 
Hamiltonian has two degrees of freedom, one first integral other than H as- 
sures the Liouville integrability of H and if the submanifolds of the foliation 
are compact, then they are diffeomorphic to tori [1]. 



2 Poincare sections 

For two-dimensional completely integrable Hamiltonian systems the phase 
space is foliated into two-dimensional manifolds and each integral curve of 
the system belongs to one and one only of the leaves of the foliation. In this 
case, the intersections of any integral curve with a given plane of the phase 
space ("Poincare sections" or "phase sections") consist of points arranged 
into curves, determined by the intersection of the two-dimensional Liouville 
torus with the section plane. Naturally, the Liouville tori are only diffeo- 
morphic to the standard torus, then they appear twisted and folded into the 
phase space so that their intersections with a section plane can produce sev- 
eral disconnetted closed curves. This behaviour is a characteristic feature of 
two-dimensional integrable systems. Poincare sections of non-integrable sys- 
tems, instead, can show both points arranged in curves or points shattered 
into the section plane, depending on their different initial data (Figure [T|, 
jl8j . If one want to show the integrability of a system by means of Poincare 
sections only, it would be necessary to produce sections for all possible sets 
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of initial conditions all showing points forming curves; however, this is clarly 
impossible in practice. The only proof of the existence of first integrals, and 
complete integrability, is the knowledge of their mathematical expression it- 
self and the validity of Poisson equations. Such a knowledge can be extremely 
hard to obtain for Hamiltonian with not simple potentials. An experimental 
approach is however possible: numerical integration of the integral curves of 
the system can be performed for generic initial data and Poincare sections 
can be obtained. Several computer algebra systems can do the task and 
we used the "poincare" procedure implemented in Maple 9.5. By analyzing 
the Poincare sections so obtained we can have an experimental evidence of 
integrability, in the case all sections we obtain for generic initial data show 
the features expected in case of integrable systems, and therefore focalize the 
efforts in search of a rigorous proof. We produce below some of these exper- 
imental evidences for three potentials with the symmetries of the Platonic 
polyhedra. 



3 The Calogero system and dihedral symme- 
try 

An integrable system of great importance in mathematical physics is the 
Calogero-Sutherland-Moser system [3], which is studied both in its classical 
and quantum form. This system, in its simplest form, describes the reciprocal 
interactions of three points on a line, denoted by their positions x* with 
respect to some origin, with potential 

1 1 1 



2 ■ 



It is possible to write the Calogero system as a one-point system in with 
cylindrical coordinates (r, '0,^) (see Section 8 and [5]) which, because of the 
conservation of the momentum, reduces to a two-dimensional system in the 
plane (r, with potential 

(2) V ^ 



r^sin S?/' 



It is now evident the dihedral symmetry of the system: the plane is divided 
into six identical sectors where the dynamics is the same, and the symmetries 
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of the potential are those of the hexagon (the discrete rotational symmetries 
of the regular polygons and their refiectional symmetries are called dihe- 
dral symmetries). Particles moving on the plane under the potential V are 
trapped into each sector by the infinite value of the potential attained on 
the boundaries of the sector. In [3] is conjectured the super-integrability of 
systems in the plane with potential 



sin^ kip 

with k integer and, for odd k, the expression of the corresponding constant 
of the motion is given (the proof will be published soon), moreover, these 
systems are shown to correspond to three-body interactions among three 
points on a line in the same way of the Calogero system. For these systems 
the dihedral symmetry corresponds to the symmetry of the regular polygon 
with 2k sides. Other studies about systems with dihedral symmetry recently 
appeared are [25 , [25] , [16] . The hexagonal symmetry of the potential ^ is 
directly related to a cubic in the momenta first integral of the system and 
the same holds for more general potentials [S]. It must be remarked that 
the hexagonal symmetry is completely hidden if the potential is represented 
as ([T]). The main idea of the present investigation is (try) to show that, in 
analogy with the Calogero system, finite symmetries of the potential, maybe 
together with some other ingredient, lead to first integrals. 



4 Systems with Platonic symmetries 

It is natural try to generalize the previous results to the three dimensional 
space, then, to search for integrable Hamiltonian systems with potentials of 
the form 

(3) V ^ 



p'fiO,iP) 



in spherical coordinates (p, 6, ip) where / is a function on the sphere , 
with polyhedral symmetry. Integrable systems of this kind will be super- 
integrable if embedded in (with 5 independent constants of the motion) 
and equivalent to superintegrable interactions among four points on a line, 
as we show in Section 9. The resulting natural Hamiltonian system in the 
three-dimensional Euclidean space is sometimes referred to as "conformal 
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mechanics" and its construction is generalizable to the building of integrable 
and superintegrable systems from spheres to n-dimensional Euchdean 
spaces [5] , [6] , [10] (where connections are made with Calogero system, Higgs 
oscillators and supersymmetric mechanics). 

We limit our investigation to potentials with the same symmetries of the 
five Platonic polyhedra as the simplest three-dimensional generalization of 
the dihedral symmetry manifested by the Calogero and the other systems 
seen above. There is no need to remark the relevance of Platonic polyhedra 
in philosophy, arts and science, more can be found for example in [7]. The 
Platonic polyhedra are those polyhedra whose faces are all made by the same 
regular polygons and whose vertices belong to a sphere. Their symmetries 
are the rotations and reflections leaving fixed the center of the polyhedron 
and making the faces to correspond, leaving in this way unaltered the ap- 
pearence of the solid: the polyhedral symmetries are those rotations and 
reflections which leave invariant the polyhedron. Since Euclid's times it is 
known that there are only five possible Platonic polyhedra: tetrahedron, ex- 
ahedron (or cube), octahedron, dodecahedron and icosahedron. Actually, 
esahedron-octahedron and dodecahedron-icosahedron form so called dual (or 
reciprocal) couples (the centers of the faces of one correspond to the vertices 
of the other) and each couple share the same symmetries while the tetra- 
hedron is self-dual. The same symmetries are shared by the Archimedean 
polyhedra obtained from the Platonic ones, and by their duals the Catalan 
solids [7], [1], [23]. Namely, the symmetries we are considering are the poly- 
hedral groups denoted by T12, O24 and Iqq (isomorphic to ^4(4), 5(4), A{5) 
respectively, A{n) and S{n) denoting the alternate and symmetric grups of 
degree n respectively), where the lower number is the order of the group. 

An example of superintegrable system with platonic symmetry whose 
first integrals can be explicited is given in [TT] . It is obtained from the three 
particle Calogero system D3 [20j and interpreted as potential generated by 
six centers of force on the sphere, each one on the vertices of a cuboctahedron 
(one of the Archimedean polyhedra) . It is a superintegrable system with the 
symmetries of the exahedron-octahedron, and its potential on the sphere is 

9(8-tan^g)^ 12 

~ 2(3 tan^ 6-8 + tan^ 6 cos Stpy ^ 3 tan^ 6-8 + tan^ 6 cos 3V^ ^ 

9 

^4sin2^(l + cos6^)' 
its first integrals can be obtained by following [TT]. Here, the integrability 
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of the system is derived from that of the original Calogero system D^, the 
complexity of the potential make almost impossible any direct inquiry in 
that direction. We try here to build simpler integrable systems with platonic 
symmetries in order to provide more suitable matter for further analysis. 

Symmetry groups for Platonic polyhedra and their polynomial invariants, 
i.e. polynomials in CArtesian coordinates [x, y, z) left invariant by the sym- 
metry groups of the corresponding polyhedra, are well known. Examples of 
applications of platonic symmetries in physics and mathematics are [2] or 
[23l| . generalizations are currently object of investigation relatively to more 
general symmetry groups |T2]. 

It can be shown that all invariant polynomials over the reals M (or the 
complex C) for each one of the symmetry group of the Platonic polyhe- 
dra can be written as polynomials, over the field M (or C), in the variables 
[[/i, U2, t/3], where the Ui are suitable homogeneous poynomials in Cartesian 
coordinates {x,y,z) (an instance of the celebrated Hilbert finite base theo- 
rem) [5], [2T]. Because of their homogeneity, these polynomials, written in 
spherical coordinates, can always be factorized into the form Ui = p''fi{9,ip) 
for some positive integer k, this correspond to some "conformal invariance" 
of the system [lOj. Therefore, the functions fi on so determined carry all 
the polyhedral symmetries of the original polynomial, the same do f^^ and 
any function of the form g{p)f^^ for arbitrary functions g{p)- In analogy 
with the two-dimensional dihedral case, we will consider potential functions 
on §^ of the form V = and in of the form W = p~'^f~^. To pro- 
duce our examples we consider some of the polynomials of the bases (f/j) as 
given in p5], due to the fact that they have possibly the simplest algebraic 
expression for functions with the desired symmetries. 

In these bases, the lowest-order polynomial is always x"^ + y"^ + which 
obviously is invariant under all rotations and reflections leaving fixed the 
origin, while, 

T = xyz 

for the tetrahedron, 

/o 222 
= xyz 

for the exahedron-octahedron and 

/ = -z{2x + z){x^-x'^z^ + z^ + 2{x^z-xz^)+5{y'^-y^z'^) + 10{xy'^z-x^y^)) 

for the dodecahedron-icosahedron, are characteristic of each of them. The 
polynomial TO = x'^y'^ +x'^ z"^ +y'^ z^ is common to the bases of tetrahedral and 
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octahedral invariant polynomials. By writing these polynomials in spherical 
coordinates and by factorizing out the radial terms, we obtain, as described 
above, the functions fi for the polynomials T, O and / and from ff^ the 
following "platonic" potentials on the sphere: 

Vt = (sin^ 9 cos 9 cos sin tp)~^ 

with tetrahedral symmetry, 

Vo = vS 

with cubic-octahedral symmetry, and 

Vi = — cos~^ 9[cos^ 9 — 5 sin^ 9 cos^ 6^ + 5 sin"^ 9 cos 9+ 
sin^ ^(32 COS?/' sin^ ip — 24 cos ■?/' sin^ ip + 2 cosip)]~^, 

with dodecahedral-icosahedral symmetry. 

Still, it is not easy to find directly the expressions of first integrals of 
these systems, even by using computer algebra methods as we did in the 
lower-dimensional case Therefore, we integrate numerically the natural 
Hamiltonian systems with the three potentials of above on the sphere and 
analyze their Poincare sections. 

Poincare sections are meaningful only if orbits belong to a compact sub- 
manifold of the phase space, in this case each orbit, numerically computed, 
winds itself in general several times around the Liouville torus and can pro- 
duce several points on a suitably chosen section plane. This happens in our 
case in the regions of §^ where the potential is positive. In the remaining 
regions the potential generates a force pulling the particle towards the bor- 
ders of the region and the orbit crosses the section planes only few times. 
The structure of the Hamiltonian in the negative-potential regions could be 
studied by considering — V instead of V, which of course are different systems. 

It is remarkable that in all these examples the sphere is partitioned by the 
lines of the zeroes of the fi into separated regions. Because on these lines the 
value of V goes to infinity, particles moving under the potential V cannot 
cross the borders of these regions. In these regions, the potential admits 
just one critical point which is either a maximum or a minimum (Figure |2|, 
therefore, the structure of the dynamics is essentially simple: there are on the 
sphere dynamically separated regions each one with some simple dynamical 
structure. The same happens for the Calogero system considered above: the 
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sphere is partitioned into six regions by the infinities of the potential, 
regions equivalent under the symmetries of the hexagon. 

Possibly as a consequence of that additional structure, the corresponding 
Hamiltonian systems appear to be Liouville integrable: the computation 
of Poincare sections for several randomly chosen initial conditions in the 
regions where the potential is positive shows in fact always the curvilinear 
features of the intersections of the integral curves with the section planes, 
giving evidence of an independent first integral K{p0,p^,9,ip) at least for 
each system (Figure |3|. This happens also for the systems with potential 
—V in the regions where V is negative. The simple potential in each region 
can allow the existence of a local first integral which is extended on the 
other regions corresponding under the platonic symmetry group. The same 
approach with potentials on the sphere obtained from the fi instead of the 
fj^^ of above do not show yet signs of integrabihty. Some computations made 
with the potential obtained in a similar way from the third polynomial of 
the icosahedral base (Ui) given in [5] indicate its integrabihty, however, much 
less sections have been considered in this case. 

Not all the potentials with platonic symmetry obtained from the polyno- 
mials Ui seem to be integrable; as an example consider 

Vto = sin~^ 6'[cos^ ip — cos'' 'ip — cos^ ip cos^ 9 + cos"^ ip cos^ 6 + cos^ 

obtained from the polynomial TO of above, which shows six isolated points of 
maximum with infinite value there. In this case the sphere is not partitioned 
into regions with a simple behaviourhed potential. Neither all integrable 
potentials with platonic symmetries partition the sphere into dynamically 
separated regions with a simple potential within, for example, Vco- 

The integrable dynamics just analyzed on the sphere can be extended 
to the three-dimensional space with the potentials W = p~'^V . The new 
natural Hamiltonian systems are integrable in if the original ones on the 
sphere are, and their embedding in dimension 4 is superintegrable as we 
prove in Sections 7 and 8. An harmonic term proportional to can be 
added to each one of the previous potentials in keeping the integrabihty 
and allowing for finite trajectories of the systems. Then, integrabihty can 
be studied by analyzing the orbit structure of the system by following the 
approach developed in dimension 2 by [23]. Even if we do not know yet the 
expression of the possible first integrals, likely they are polynomial in the 
momenta. 
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In Sections 8 and 9 we see as Integrable systems on lead to super- 
integrable 4-points systems on a line, in the same way integrable systems 
on S"~^ can be interpreted as superintegrable n-body systems on a line [5]. 
In this perspective could be interesting to analyze higher dimensional poly- 
topes (by using a Coxeter's expression, a polytope is the general term of 
the sequence: point, segment, polygon, polyhedron,...), in Euclidean or non 
Euclidean spaces P]. 

5 Remarks 

The partition of the configuration manifold into "simple" dynamically sep- 
arated regions united to some suitably rich symmetry group on the same 
manifold seems a good indicator of integrable potentials. For example, on 
the Euclidean plane potentials of the form 

Vi = (sin"/ixcos''%)~\ 

with a, b, h, k positive integers and (x, y) cartesian coordinates, partition by 
their lines of infinities the plane into rectangular regions where the potential 
has only one point of maximum or minimum and a simple behaviour. In each 
region of the partition the Poincare sections show a dynamics compatible with 
integrability. On the plane, these systems admit translational symmetries 
along and y of magnitude ^vr (resp ^tt), where Sx (sy) are even 

integers for a (resp. b) odd and any integer otherwise, all generated by 
a I (^ for odd a) translation in the x direction and by a | (^ for odd 
b) translation in the y direction. The lines of reflectional symmetry have 
equations x = and y = ^2F^^- seems likely that integrability holds 
also for the three-dimensional natural Hamiltonian systems in the Euclidean 
space with potential of type 

Wi = {sin'^ hx sin'' ky sin'^ I z)~^ 

with c and / positive integers, as well as for the generalization of these systems 
to higher dimensions, where again the space is honeycombed into dynamically 
separated cells with simple potentials in each of them and the symmetry 
group is generated by the elementary translations of magnitude f , f , f (twice 
in the cases of odd a, b, c respectively) in the three coordinate directions, 
plus the obvious reflectional symmetries. In all these cases, assuming they 
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are integrable, the integrability in each cell, which is due to some local first 
integral, is extended to the whole configuration manifold by the symmetry 
group on it and the local first integrals become a global one. 
Potentials of the form 

V2 = {sin" he sin'' ktjj)-^ 

instead, have polyhedral but not platonic symmetry and partition the sphere 
into a meridian-parallel web of "simple" dynamically separated regions by 
their lines of infinite value. In this case, the rotations of the dihedral symme- 
try group on the configuration manifold are generated by the rotation around 
the poles of amplitude | for b even and | for odd b, a dihedral group. The 
planes of refiectional symmetry are defined by the integer multiples oiijj = 
for b even and = 2| for odd b. There is the additional refiectional symme- 
try with respect to the equatorial plane when a is even or h is odd. While 
appearing in general integrable, in some cases their Poincare sections seem 
to show integrable behaviour only in the regions which are not neighbouring 
to the poles, for example for a = 4, h = 2, b = 1, k = 3. 

A behaviour analogous to the systems with potentials of type V2 is shown 
in the Euclidean plane by the potentials of type 

V3 = {sin'' hr cos'' k^)-^ 

with, mutatis mutandis, all possible permutations of sin and cos functions, 
where (r, 6) are polar coordinates. In the regions neighbouring to the origin 
they seem do not to be always integrable, differently from the others regions, 
for example when a = b = 2,h = k = l. Here, an evident dihedral symmetry 
group exists. In all these cases, some kind of local first integral could exist on 
the quadrangular regions only, but not on the triangular ones neighbouring 
to the origin. A local integral is extended by the symmetries of the config- 
uration manifold to a semiglobal first integral on the whole manifold minus 
the regions surrounding the poles of the coordinates: the whole plane minus 
a disc centered on the origin of the coordinates. Due to the dynamical sep- 
aration of each region, systems of this kind could be considered completely 
integrable as far as initial conditions are not chosen to be into the "bad" 
regions around the poles. 

Differently from these last examples, the Calogero system considered 
above, the systems in [5] with dihedral symmetry and the particular case 
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of V2 with a = 2, /i = 1 on the sphere, admit the additional symmetry of sep- 
arabihty. Indeed, in all these cases the differential equations of the dynamics 
are separated for the two coordinates thanks to a quadratic first integral (see 
|5] for the first two cases). On the sphere, the natural Hamiltonian with 
potential 



sin^^' 

^2 



for any function F(ijj), admits the quadratic first integral + 2F{^) and 
is therefore completely integrable and separable as in the previous cases. 
Symmetries like the latter are sometimes called "hidden", because not im- 
mediately recognizable from the expression of the potential, as instead is for 
example the central symmetry for the Kepler system which is associated with 
the conservation of the angular momentum. In some cases, hidden symme- 
tries can be unveiled by much more evident ones. It is the case of the three 
body Calogero system shown above where the hidden symmetry is a third- 
order polynomial in the momenta, unveiled by the hexagonal symmetry of 
the potential written in the center of mass frame. Remarkably, in [15] the 
three-dimensional system with Hamiltonian 

H = p2 I I pI I « I 1 ( Pi + ^ 

^ p^siv?h9 p ycos"^ h6 siv? hO cos^ kif) sm^ hO siv? kip 

where a, /3j are real parameters, is shown to be maximally superintegrable for 
all /i, k rationals. If reduced to the sub manifolds p = const., for a = (3i = 
and h = 1 the system becomes an instance of the natural Hamiltonian on 
the sphere with potential V4 and dihedral symmetry. 

We think to have provided here several evidences for the presence of 
hidden symmetries connected both with platonic symmetries of the potential 
and with the partition of the configuration manifold into simple dynamically 
separated regions. Anyway, all the examples of above are just hints to further 
inquiries and we do not pretend to give here a detailed description of their 
behaviour. 



6 Quantization 

The classical systems here considered can be transformed into quantum me- 
chanical ones by standard quantization techniques. Indeed, the Calogero 
system is in origin quantistic. Instead of Hamiltonians we have in this case 
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Schrodinger operators on the sphere or on the Euchdean three-dimensional 
space, instead of Poisson commuting quadratic in the momenta first inte- 
grals, commuting second order differential operators. More subtle is the 
quantization of higher-order first integrals [8j or the definition of quantum 
superintegrability • We do not consider further here the quantum version 
of our systems. 

7 First integrals in 

The procedure described below applies to every Hamiltonian integrable sys- 
tem on the sphere with potential of the form ([s]), in particular to the systems 
described in Section 4, provided they admit a first integral and are, conse- 
quently, completely integrable. The natural Hamiltonian on §^ is 

2 sm U ^ 

and our platonic systems correspond to the case when V is either Vr, VooiVj. 
Let us call H2{O,ilj,p0,p^) an independent first integral of Hi, (the unknown 
first integral inferred from the structure of the Poincare sections in our case). 
The system is then Liouville integrable. Let, in spherical coordinates (p, 6, ip) 

of M^ 




the natural Hamiltonian with potential W = + ^p^ where the original 
potentials on the sphere are modified by a harmonic term with parameter 
k G M^. The system determined by H3 is completely integrable if Hi is. 
Indeed, the three functions H3, Hi and H2 are functionally independent and 
all in involution if Hi and H2 are. The equations of the motion are partially 
separated into equations in Oyip, for the dynamics projected onto the sphere 
of fixed radius po (the dynamics of ifi), and 

2 

pI = 2/i3 - kp^ -h2 

where hi are the values taken by Hi in the given initial conditions. 
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8 First integrals in M^. 



The Hamiltonian of above can be extended to IR4 [S], [in] by introducing 
cylindrical coordinates (m, p, 9, and obtaining the natural Hamiltonian 

H, = ^ipl + 2Hs) 
which admits the trivial first integral 



and, if = 0, the less trivial 

1 

A p 

in these coordinates, p is the distance from the axis of the three-dimensional 
cylindrical hypersurfaces. The functions Hi, H2, H4 and are still inde- 
pendent and in involution with each-other, while Hq, if k = 0, is in involution 
with Hi, H2, H4 and independent from Hi, H2, H^ and H^. Therefore, the 
system of Hamiltonian H^ is (minimally) superintegrable in M^. Other ex- 
amples of this kind of extension, apphed to superintegrable Evans systems, 
are given in [6j. 



9 Four points on a line 

The system of above can be interpreted as a natural Hamiltonian system 
describing reciprocal four-body interactions among four points on a line by 
the change of variables 

= , ^ (x^ + ...+x^ - 7 = 1... 3 



u 



4 



1 



where = x, = y, = z, = u and denote the positions of four 
points on a line whose dynamics is still described by the Hamiltonian 
and keep the same integrals of the motion. The integral H^ correspond to 



14 



the conservation of the momentum of the system of four bodies. Evidently, 
the previous identification between one-point systems in and four body 
systems on M is essentially unaffected by phase shifts in 6 and ■j/^, namely 
transformations 6 ^ 6 + 6q and iIj ^ ip + ifjQ, because such phase shifts do 
not modify the dynamics of the four points but only their reciprocal position 
on the line. An example of the equivalence just remarked is given in [5] 
between the Calogero and Wolfes potentials. The procedure can be extended 
to n points on a line [6] . 

10 Numerical procedures 

The Poincare sections of this paper have been obtained by using the "poincare" 
procedure implemented in Maple 9.5, which is a fourth-order Runge-Kutta 
algorithm. For each section the integral curves have been numerically inte- 
grated for a Hamiltonian parameter t ranging typically between ±100 or ±50 
and several hundreds of crossing points are obtained. The section planes par- 
allel to or (pi,p2) have been selected in order to show intersections 
for all the four sets of initial conditions and the choice is not particularly 
relevant regarding the shape of the sections. The procedure "poincare" au- 
tomatically check the maximal percentage deviation from the initial value 
of the Hamiltonian along the computed integral curve, giving in this way 
some measure of the confidence of the integration. If the percentage devi- 
ation exceed some given maximum, the accuracy of the integration can be 
increased, for example, by refining the discretization of the t interval. In our 
computations we allowed deviations up to 1 x 10~^% of H, even if in most of 
the graphics the deviation is typically 1000 times smaller. Point crossing the 
section plane are assumed do not describe curves if, after further refinements 
of the t discretization, they keep that behaviour. Obviously, the curves in the 
section planes are obtained after some minimal discretization and are stable 
for finer discretizations. The practical minimum of the discretization step on 
our computer for t in the intervals of above is around 0.002. 

11 Conclusions. 

In this paper we produce evidences of complete integrability for two dimen- 
sional Hamiltonian systems on the sphere with the symmetries of the pla- 
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tonic polyhedra. By assuming the effective existence of a first integral, we 
sliow llow to extend tlie systems from tlie splierc to an intcgrablc system 
in tlie three-dimensional space and how to build superintegrable systems in 
dimension four corresponding to four-body interactions on a line. Even if 
no explicit first integral of the platonic systems on the sphere is obtained, 
its existence seems more than probable. The integr ability of these systems 
seems ascribablc to the presence of a partition of the configuration manifold 
into dynamically separated regions, each one with a simple structure of the 
potential allowing a local first integral, and of a symmetry group allowing 
the extension of the local first integrals into a global one on the whole con- 
figuration manifold. The conjecture is extended to a class of similar natural 
Hamiltonians in Euclidean spaces. If the extension fails in some of the re- 
gions, the system can be considered integrable in a scmiglobal sense. The 
approach seems fertile for further inquiries: first, towards a determination of 
the first intergrals, second, towards the extension of the approach to higher 
dimensional polytopes and to non-euclidean spaces. 
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black: H = 3.7BI245000; red: H = 4,151683000; 




Figure 1: A cross section showing intersections of two distinct trajectories for Vro- 
Black squares correspond to an orbit lying on some 2- dimensional suhmanifold of the 
phase space while the magenta ones correspond to an orbit not lying on a surface. 
Therefore, the system do not show behaviour of complete integrability. 



19 



Figure 2: For the potentials Vt, Vq and Vj with tetrahedral, cubic- octahedral and 
dodecahedral-icosahedral symmetry respectively, some isopotential lines are drawn on 
the sphere Aquamarine- green-blue denote decreasing negative and magenta- orange- 
red increasing positive values of the potential. Black lines denote infinite values of the 
potential, therefore, they determine regions of the sphere where the motion is confined. 
Regions with identical behaviour of the potential correspond under the symmetry group 
of the associated polyhedron. The isopotential lines show the simple structure of the 
potential in each region of the sphere bounded by the black lines. 
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Figure 3: Are shown here examples of Poincare sections, with qi = 6, pi = pg the 
momentum conjugate to 6, of integral curves of natural systems on the sphere with 
potential, from left to right, Vr, Vq and Vj with four distinct initial condition sets, dis- 
tinct by their colors. After some minimal level of accuracy in the numerical integration 
of the integral curves, their intersection points with a plane {pi,q^) in the phase space 
describe closed curves, a behaviour compatible with complete integrability. 
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Figure 4: Poincare sections in planes parallel to (pi = pe,P2 = Pi/') /^'^ ^o, Vi 
respectively, with the same initial conditions of Figure \^ 
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